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Abstract 

We propose that the development of odd-frequency superconductivity is 
driven by the formation of neutral three body bound-states. Using a three- 
body bound-state ansatz we develop a mean-field theory for odd-frequency 
pairing within the Kondo Lattice model. Three body bound-state formation 
leads to the formation of a gapless band of fermions with a neutral, spinless 
Fermi surface. We discuss the low energy excitations of these modes, sug- 
gesting them as a possible explanation for the absence of anisotropy in the 
thermal conductivity of heavy fermion superconductors. 
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Although three body fermion states are ubiquitous in many branches of particle, nuclear 
and atomic physics, there has been very little study of their impact on collective condensed 
matter behavior. In this paper, we propose that three body bound-state instabilitieJil can 
precipitate a phase transition, and we identify the broken symmetry state that forms as 
an odd-frequency superconductor. Berezinskiii originally proposed that superconductivity 
could occur through the development of a state with an odd frequency gap function 



Unlike conventional superconductivity, there is no weak-coupling Cooper instability asso- 
ciated with this type of state.il Indeed the proper identification of an appropriate class of 
many body instability for odd-frequency pairing has proved a major obstacle to concrete 
theoretical realization of Berezinskii's proposal. 

Our particular interest in odd-frequency pairing is motivated by heavy fermion 
superconductors,&i where superconductivity appears to involve the active participation of 
local moments in the condensation process. Vanishing coherence factors particular to odd 
frequency pairing may offer a means to reconcile the large linear specific heat and isotropic 
thermal conductivity observed in these superconductors with the absence of a linear Kor- 
ringa relaxation in the NMR data.iHl2l In this paper we show how the development of such 
hypothetical states can be directly associated with the formation of three-body bound-states 
between electrons and local moments. 

A three-fermion bound-state is a composite fermion where selection rules stabilize it 
against decay into its constituents. A good example is the IIe-3 atom: a neutral bound- 
state formed from a 5 = 1/2 nucleus and two electrons. Since the low-energy dynamics of 
the atom involve rigid center-of-mass motion, correlation functions of the bound fermions 
are determined by their factorization into triplet contractions 



where ^"^{x) creates the He-3 fermion, i/) represents the electron fields. A/" the nucleus, and 
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A is the atomic wave function of the He-3 atom. It is this factorization into "three body 
contractions" which enables us to replace the composite operators by a single fermion $(a;) 
in the low energy physics of He-3. 

Can such bound-states form collectively in an electronic system ? A well-known theorem 
of Yangini precludes the development of any literal off-diagonal long-range order of a three- 
body fermion. Here we propose an alternative possibility, where two electrons and a hole 
bind to form a neutral fermion through the collective formation of a three-body condensate:0 



1 \ 1 r 

iJJ{l)^Pp{2)^P^{3)= j Kp^{l,2,?>-x)<P{x)dx (3) 

Here = (p^x) represents the neutral bound-state fermion. Aq,/?^ is a complex wave func- 
tion which carries the phase information associated with the charge of the bound electrons, 
playing the role of a collective variable. We shall argue that if the phase of the three-body 
wave function develops a rigidity, an odd frequency superconductor is formed. For simplicity 
we consider a wave function that is symmetric in electron co-ordinates 2 and 3, and hence 
antisymmetric in the corresponding spin variables 

A,^^(l, 2, 3; x) = aiy^^A^il, 2, 3; x) (4) 

where cx^ is a Pauli matrix. The coarse-grained average A{x) = 2 3 ^(l? 2, 3; x) plays the 
role of a spinorial order parameter.^! By contracting the spin-indices, we see that A describes 
the binding between spins and electrons. 



[3(1) ■ a^p]M'2)= / dxA^il, 1, 2; (5) 



where S = ^ip^aip is the spin density. 

A specific model for this type of bound-state instability is provided the "Kondo lattice" 
Hamiltonian 

^ = Eek<^^k + E^»*b1 (6) 
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where ip^^ is a conduction electron spinor, coupled to an array of 5" = 1/2 local moments Sj 
via an antiferromagnetic exchange interaction 

Hintlj] = Ai^^j^i^j) ■ Sj. (7) 

where ipj denotes the conduction electron in a tight binding representation. This Hamilto- 
nian is of particular interest as a toy model for heavy fcrmion metals. An electron scattering 
off a magnetic ion at site j couples directly to the three-body spinor 

0« = i^j ■ ^ap)i>jl3 (8) 

Irreducible Green functions of ^ determine the self-energies of the conduction electrons. In a 
superconducting state, the normal and anomalous components of the electron self-energies 
are given by 

where double brackets denote one-particle irreducibility with respect to conduction electrons 

— * 

and K = {ujjk). The development of bound-state resonances associated with ^ can thus 
severely modify the electronic properties of the lattice. 

To provide a controlled treatment of fermionic bound-states in the Kondo lattice model, 
we develop a,l/N expansion within a class of "0{N) Kondo models" . In this generahzation, 
the local anticommuting algebra of the original spin-1/2 operators is preserved 



{Sj, Sj} — 5^1/ 



1 (10) 

but the number of components of each operator is generalized from three in the model 
of interest, to an arbitrary odd number N. This is the algebra of N component gamma 
matrices. For general N, the conduction electron spinors have dimension 21^/^1, where [ ] 
denotes the integer part. For all N, the Kondo interaction Hamiltonian can be written 

H^nt\J] = -^{^jQ (11) 



where we replace J = J/N. 

We now calculate the correlation function of ^ in the normal state, in a scheme that is 
exact to order 0{1/N). Since ^ is a fermion, a fluctuating component of the operator must 
survive development of bound-state singularities in its correlation functions. To describe 
these intrinsic fluctuations we introduce the device of a "spectator fermion" {(pj): a neutral 
fermion (pj = (pj^ defined at each site which which satisfies canonical anticommutation 
rules, = 5jk, commuting with the spin operators and anticommuting with all other 

fermions. 0^ , so defined, commutes with the Hamiltonian and decouples from the dynamics; 
its correlation functions display intrinsic fluctuations 



(T0,(l)0fc(2)) = ^sgn(ri-r2) (12) 

We construct a "bosonic partner" to the three-body operator ^ as follows 

i, = 20,^,. (13) 

Since 0| = |, we may invert this relation to obtain S,j = (pj^j- The spectator fermion factors 
out of the dynamics, so we may relate the correlation functions of Aj and ^j, 

(i,(r)it,(0)) = 2sgn(r)5,,((e,(r)e/(0))) (M) 

We can also reversibly transform the spin operators into fermionic counterparts 

Vj = '^<i>jSj (15) 

which satisfy {^^'^[j], ?7^[A;]} = S"'^6jk and are realrj°-[j] = //"^[j]- The canonical commutation 
properties of these "Majorana" spin fermions permits us to use diagrammatic methods 
without invoking a Gutzwiller projection.0 Since Aj = [ffj ■ (T]ipj, the original interaction 
can be written 

H^nt[J] = -^A^^A, (16) 

Consider the diagrammatic expansion for the local susceptibility 



(A(r)At,(0)) = 5,,x(r) (17) 

in powers of J /N. (Fig. (1.)) The non-interacting value is 'X^{t) — ySgn(T)^(T) where 
Q{t) = (^j(r)^''"j(0)) is the local electron propagator. In the frequency domain 

tanl^V?) ,^__^2.T„) (18) 

is logarithmically divergent at low energies x°(0) ^ Npln[D/T], where p is the conduction 
electron density of states and D the bandwidth. Each fermion loop contributes a factor N 
and each interaction hue contributes a factor 1/N, thus the leading order 0{N) contribution 
to X is given by the RPA sum 

XRPA{^^n)-N (19) 
(1 - txO(^i^n)) 




Fig. 1. Diagrammatic expansion of x(iz^n) in a expansion for the 0{N) Kondo 
model: solid lines refer to conduction electron propagators, dashed lines to spin fermions, 
and a wavy line denotes the interaction vertex. 
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where Nx'^ = If we now transform back to the corresponding fermionic response function 
of 1^, we find 

UrWm = lsgn{T)xRPA{r) (20) 

Since x°(0) is logarithmically divergent, this function develops a singularity at a finite tem- 
perature Tc ~ De"^/'''^, where 1 = ■^x°(0). This singularity signals the formation of a three 
body bound-state. 

We now apply our bound-state ansatz to the three body spinor ^j, writing 

e,(t) = (A,)0,(t) + 5e,(t) (21) 

where the first term is the three body contraction of the operator. This ansatz implies 
long-time correlations of the three body operator 

((0(^)e%(0))) ,— - %sgn(r)A, ® A/ 

((0(^)6(0))) ^sgn(r)A, ® 



For = 3, The Fourier transform of (|2^), inserted into (|^) then gives 

= ® V^mK)<P{-K)) = IIK^Z!!. (23) 

UJ 

where Vj = ^{Aj). This establishes the direct link between three body bound-state forma- 
tion and an odd frequency pairing amplitude that scales as 

To develop the mean-field theory for the condensed phase, we substitute the bound-state 
ansatz into the Kondo exchange interaction, so that Hint[j] — > Hint[Vj], where 

H,nt[V,] = 2[(e^-0,r, + (H.C.)) + + 0(5^^50 



Using dll), 2^t0 = ^t(^ . thus 



Hint[Vj] 



Y!^ (24) 

Remarkably, the fusion of the neutral bound-states with the the spin variables transforms 
them into propagating spin fermions. This form of mean-field theory was previously derived 
using an abstract Majorana spin-representation of spin-1/2 operators.! 



Formation of a gapless band of neutral fermions follows naturally from our original ansatz. 
If we write Vj = -^Zj, where zj is a unit spinor, and decompose ipj into four Majorana 
components ipj = -^[i^j + ii^j ■ then the scalar component ipo decouples from the the 
hybridization term Hint[Vj] = —iVipj ■ fjj. The vector components of the conduction sea 
hybridize with the spin fermions to form a gap ~ V'^/D. Residual "scalar" components 
of the conduction sea thus split off below the quasiparticle continuum to produce a gapless 
band. (Fig. 2.) Spin and charge operators are off-diagonal bilinears in this Majorana basis, 
which ensures the neutrality of the gapless Fermi surface and leads to coherence factors that 
vanish linearly with energy.i 




Fig. 2. Schematic illustration of the excitation spectrum predicted by the mean-field theory. 
Formation of the three body bound-state splits a gapless band of neutral fermions off from 
the charged quasiparticle continuum. 

The use of spectator fields differs fundamentally from the more familiar use of "slave" 
fields, since the combined Hilbert space of particles and spectators replicates the physical 
Hilbert space. There is no unwanted subspace with unphysical matrix elements to be pro- 
jected out. In a path integral language the decoupling of the spectator fields is associated 
with a unique fermionic gauge invariance. The appropriate Lagrangian is 

C = Y^i^.^dr^, + ^ff.dff, + Hi,,t[Vj]} + (25) 
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where Vj is now a fluctuating bosonic spinor field. This Lagrangian is derived from the 
observation that the combined Hilbert space of spectator fermions and spins is equivalent to 
the Hilbert space of the Majorana spins: Tr^ = ^r^^. The saddle point approximation to 
the path integral of this Lagrangian corresponds to the mean-field theory described above. 
Suppose we couple a Grassmann source term aj to the (p field 

C + Y,aj(pj, {aj = draj ) . (26) 

j 

For N = 3, (pj = —2irjjrj^rjjM We can always gauge away such a term by a transformation 

Vj ^Vj + 2ajSj (27) 

under which the spins Sj and the Hamiltonian H = J2j Hint[j] are invariant. This is a 
"supersymmetric" transformation, generated by the (pj fields in the unitary transformation 
U = Y[j{^ +(yj(pj)- Splitting the Berry phase term into separate contributions from the spin 
{Cb[S]) and field, 

\ffdrff = Cb[S] + ^0(9^0 

. ? (28) 

/:b[5] + -(0 + «)5.(0 + «) 

we find that this term transforms as 

J dT[-f]jdrf]j + a^cPj] j dT[-f]jdrf]j - if'j^Oij] (29) 

By differentiating with respect to the source fields, we confirm that the field is decoupled, 
with a free propagator (0^(1)0^(2)) = ^ = isgn(ri — T2). 

These connections between three-body bound-state formation and odd-frequency pairing 
provide insight into its possible application to heavy fermion physics.i Coherence factors that 
are linear in excitation energy (k|p|k) oc (klS'^lk) oc uo and the development of a quadratic 
temperature or frequency dependence in a wide variety of response functions, such as the 
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normalized NMR relaxation rate ^r^, the transverse ultrasound attenuation q;^, the deple- 
tion of the superfluid density Aps and the quasiparticle conductivity (t{uj), 

(^1^, "T, Aps{T), a{u})^ oc max(T^ a;^) (30) 

are seen as a necessary consequence of neutral bound-states. These features are expected 
to co-exist with a linear specific heat and an essentially isotropic thermal conductivity. 
Crystalline anisotropy, not included in our discussion, may influence bound state formation, 
for the local bound-states should acquire a specific crystal field symmetry F of the lattice. 
The three-body spinor will take the form ^rj = • Sj)iprj where iprj = Sk7rk^*^^-''V'k 
a conduction spinor with symmetry F. Hybridization with conduction electrons will reflect 
the same symmetry — Vj^^. Nodes of the microscopic crystal field f-wave function of the 
local moment can thus lead to gap zeroes oc 7^^, suggesting an interesting possibility of 
a co-existence between neutral Fermi surfaces and gap lines of conventional quasiparticles. 

We have attempted to elucidate the physics of odd-frequency superconductivity with 
the proposal that it is driven by the formation of neutral three-body bound states. This 
hypothesis provides a physical interpretation of the appearance of Majorana spin excitations 
in the Kondo lattice model, and may be useful in developing our understanding of heavy 
fermion superconductors. 

We would like to thank E. Abrahams, P. W. Anderson, N. Andrei, L. loffe, G. Lonzarich, 
D. MacLaughlin, A. Ramirez and A. E. Ruckenstein for discussions related to this work. This 
work was supported by NSF grant DMR-93-12138, a grant from the SERC, UK and CNPq, 
Brazil. 
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Elementary defects of charge e spinors correspond to a phase change of tt, and are associ- 
ated with a flux quantum {h/e)-n = h/2e. 

For odd N, using the property Yla^"' — , we may invert the transformation, 

4>j = (— 2i)^~ Da^ab]) = (pjffj. For N — 3 this recovers the Majorana spin represen- 

— * 

tation Sj — —{i/2)fjj x ffj. 
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